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Abstract 

Our previous calculation of the spin-dependent structure function g2 is revis- 
ited. The interest in this structure function is to a great extent motivated by 
the fact that it receives contributions from twist-two as well as from twist- 
three operators already in leading order of l/Q^ thus offering the unique 
possibility of directly assessing higher-twist effects. In our former calculation 
the lattice operators were renormalized perturbatively and mixing with lower- 
dimensional operators was ignored. However, the twist-three operator which 
gives rise to the matrix element d2 mixes non-perturbatively with an operator 
of lower dimension. Taking this effect into account leads to a considerably 
smaller value of d2, which is consistent with the experimental data. 
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Typeset using REVT^ 



I. INTRODUCTION 



The nucleoli's second spin-dependent structure function g2 is of considerable phenomeno- 
logical interest. The most important theoretical tool for its analysis is the operator product 
expansion (OPE) [Q. In leading order of l/Q^, (72 receives contributions from both twist-two 
and twist-three operators. It thus offers the unique possibility of directly assessing higher- 
twist effects. The twist-three operator probes the transverse momentum distribution of the 
quarks in the nucleon, and has no simple parton model interpretation. 

In leading order of l/Q^, and for massless quarks, the moments of g2 are given by 

for even n > 2 in the flavor-nonsinglet sector. Here / runs over the light quark flavors. The 
reduced matrix elements an\jj) and dn\n), taken in a nucleon state with momentum p and 
spin vector s, are defined by 

(P' ^l^ii^i-Mnll^^' = ^^"P''^ ■■■Pi^n + traces] , (2) 

(P' *I^Hii]-Mn}IP' = ['^^'^^'^1 ~ s^lPa)p^^2 ■ ■ - Pmu + t^ccs] , (3) 

^'l^ = (^) ^7.75B,, ■ ■ ■ B,„^ - traces . (4) 

Here /i denotes the renormalization scale. The Wilson coefficients e\^, 63 „ depend on the 
ratio of scales fi'^/Q'^ and on the running coupling constant g{^J?)- The tree level values of 
the Wilson coefficients for electroproduction are given by the quark charges Q*--^-*: 

eS5 = g(^)^(l + 0(/)). (5) 

The symbol {■ • ■ } ([■■■]) indicates symmetrization (antisymmetrization) with 

The operator (|^) has twist two, whereas the operator has twist three. The twist-two 
contribution in (|l]) is also known as the Wandzura-Wilczek contribution 0. 

Note for comparison that in leading order of 1/Q^ the moments of gi are given by the 
twist -two matrix elements an '■ 

2 fdxx^grM') = \Y, e[%VQ' , gifi')) a(/)(/x) . (7) 

Both the Wilson coefficients and the operators are renormalized at the scale /i. It is 
assumed that the Wilson coefficients can be computed perturbatively. The reduced matrix 
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( f) ( f) 

elements a„ and dii , on the other hand, are non-perturbative quantities and hence a 
problem for the lattice. (Note that some authors use a different definition of a„ and dn, e.g. 
the values given in refs. have to be multiplied by 2 to agree with our conventions.) In 
the following we shall drop the flavor indices, unless they are necessary. 

A few years ago we computed the lowest non-trivial moment of g2 on the lattice . This 
calculation splits into two separate tasks. The first task is to compute the nucleon matrix 
elements of the appropriate lattice operators. This was described in detail in 0. The second 
task is to renormalize the operators. Renormalization effects are a major source of systematic 
error. An essential feature of our previous calculation was that the renormalization was 
done in perturbation theory and hence mixing with lower-dimensional operators could not 
be taken into account. In that approach the twist-three contribution turned out to be the 
dominant contribution to both the proton and the neutron structure functions. This result 
has been recently confirmed by Dolgov et al. p . 

In the meantime, it has become possible to study renormalization non-perturbatively on 
the lattice, see e.g. This approach allows us to consider mixing with lower-dimensional 
operators. If present, it will be the dominant mixing effect in the continuum limit. Since the 
twist-three operators can suffer from such mixing, we shall extend our previous work 
by employing non-perturbative renormalization. In a recent paper |^ we have started a 
non-perturbative calculation of the renormalization constants associated with the structure 
functions Fi, F2 and gi in the fiavor-nonsinglet sector. Here we consider the case of the 
structure function g2 restricting ourselves to n = 2, the lowest moment of (72 for which the 
OPE makes a statement. A preliminary version of this work based on lower statistics at a 



single value of the bare coupling has already been presented in Ref. [10 



II. RENORMALIZATION AND MIXING IN CONTINUUM PERTURBATION 

THEORY 

The renormalization of the operators which contribute to the moments of g2 has been 
studied by several authors in continuum perturbation theory |TT|. Since the more recent 



paper by Kodaira et al. [T^ (see also [jT3[) is closest to the methods applied on the lattice, 
let us briefly recapitulate the main findings of these authors. 

They consider the case n = 2 in the fiavor-nonsinglet sector and start from the operators 

i^ - 

R''/'' = - [2ij-f5-f^D^''D''^ - ij-f5-f>'D^''D''^ - ip-f^YD^^D''^] - traces , (8) 
3 

= ^9 [e^'^'^^^F^pYi^ + - traces , (9) 

i?-/^^ = im^757"L'^''7^^^ - traces , (10) 

Rl^'^ = 1 [ip^^^^D^'^YHi 9 - + V^(i ^ - m)757"D^'^7^>^] - traces . (11) 

Here Fq,^ denotes the gluon field strength tensor, which could alternatively be expressed as 
a commutator of two covariant derivatives. Due to the relation 

Rf' = -R-J^^ + R''/^ + Rl^'^ (12) 
3 
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it is possible to eliminate one of the above operators. A one-loop calculation of the quark- 
quark-gluon three-point functions with a single insertion of each of these operators reveals 
the necessity of taking one more operator into account in the process of renormalization, 
namely the gauge-variant operator 

R'^^^ = 1 [^/>757'^a^'^7^>(i p - m)^ + m)757"(9^'^7^V] - traces . (13) 

o 

Of course, in physical matrix elements neither i?eq rior i?eqi will contribute. They show up, 
however, in off-shell vertex functions and influence the renormalization factors. 

Kodaira et al. choose Ri and Rm as the physical operators. In the chiral limit m — > 
Rm is neglected, and they obtain for the scale dependence of the twist-three piece 

j\xx'gr-'-\x,Q') = (^^^J l^xx'gr-'-'ix,^,') , (14) 



where for colors and Nf flavors 



u; = , (15) 



in agreement with earlier calculations. 

Using Rp and Ri as the physical operators, one finds that in the large- A^c limit the opera- 
tor Rp dominates the renormalization group evolution of the nucleon matrix elements. This 



has been shown by Ali, Braun, and Hiller [|14|, and was rederived in the present framework 
by Sasaki [n]. 



III. RENORMALIZATION AND MIXING ON THE LATTICE 

In a lattice calculation, the first step is the analytic continuation to imaginary times, 
leading from the physical Minkowski space to Euclidean space. From now on, all expressions 
are written for the Euclidean case (for the details of our conventions see Appendix A of 
Ref. [|TB|). Hence we have to study operators of the form 

C?^,....,„=2->7.75B,,---B,>. (16) 

We shall neglect quark masses, i.e. we consider only the chiral limit. In our earlier work PJT7[] 
we have computed the renormalization constants in perturbation theory to one-loop order. 
However, it is believed that perturbation theory cannot give reliable values for the mix- 
ing with lower-dimensional operators because non-perturbative effects are expected to be 
important. 

For a multiplicatively renormalizable operator, i.e. in the absence of mixing, we can write 

OM = ZoMO{a), (17) 
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where a is the lattice spacing. The renormahzation constant Zq is fixed by a suitable 
condition. As in the continuum, we impose the (MOM-like) renormahzation condition 

tr (rR(p)rBorn(p)"^) = , tr (rBorn(p)rBorn(p)"^) (18) 

on the corresponding quark-quark vertex function in the Landau gauge. Here rBorn(p) 
denotes the Born or tree-level contribution to the vertex function. The renormalized vertex 
function rR(p) and its bare precursor T{p) are related by multiplicative renormahzation: 

TM = Z^'ZoTip) , (19) 

where = Z^{afi) is the quark wave function renormahzation constant defined as in 
Ref. i. 

As before 0, we give the nucleon a momentum in the 1-direction and choose the polar- 
ization in the 2-direction. With these choices we use the operator 

=: O^'^ (20) 

for the twist-two matrix element a2- It belongs to the representation rj^-* of the hypercubic 
group H[4:) |[T8| , p!9[| and this property protects it from mixing with operators of equal or 
lower dimension. Hence it is multiplicatively renormalizable, and the operator renormalized 
at the scale fi is written as 

0f (/i) = Z{^>(a/i)0^^>(a). (21) 
As the operator for the twist-three matrix element d2 we take 

^f2{l]4} = i (2C2{14} " ^l{24} ~ ^4{12}) 

= ^xp(j2DiD^ + 72B4B1 - l^iD2D^ - 

-\^JD{b2-\ld)2Di)i^ij 
=: C>[^] , (22) 

which belongs to the representation of H{A). The operator (|22D has dimension five and 
C-parity + and is the Euclidean counterpart of the Minkowski operator Rp- It turns out 
that there exist two more operators of dimension four and five, respectively, transforming 
identically under H{A) and having the same C-parity, with which (^) can mix: 

^1^^(^1351 - ^4354)^^ =: 0% (23) 

(71B3B1 - 71B1B3 - 74B3B4 + 74B4B3) ^ =: 0°. (24) 

We use the definition a^y = (i/2)[7^, 7^,]. 

The operator is the Euclidean analog of Ri with the field strength replaced by a 
commutator of two covariant derivatives, and O"^ corresponds to Rm- In continuum pertur- 
bation theory Rm can be neglected in the chiral limit. On the lattice, the explicit breaking 
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of chiral symmetry induced by Wilson-type fermions, which we shall use, persists even when 
the quarks are massless. For dimensional reasons, we expect that O'^ contributes with a co- 
efficient oc and hence has to be kept. The operator O^, on the other hand, being of the 
same dimension as 0^^\ mixes with a coefficient of order g"^, which should be small. There- 
fore we discard as well as possible lattice counterparts of -Req and -Reqi, which are also of 
dimension five and hence are also multiplied by a factor of order g"^. The above mentioned 
observation P^JT3| that Rp dominates over Ri in the renormalization group evolution as 
Nc ^ oo may be taken as another indication that neglecting is not unreasonable. How- 
ever, this dominance holds only in physical matrix elements and does not apply to the mixing 
with the operators -Req and -Reqi- 

So we make the following ansatz for O^^^ renormalized at the scale /i: 

= Z[5l(a/i)C[5](a) + -Z^{ai2)0''{a) . (25) 

The renormalization constant Z^^^ and the mixing coefficient are determined from the 
conditions 

tr (r^S (p)rgL(p)^) , tr (rgL(p)rgL(p)t) , (26) 
tr (rgi(p)r5_(p)t) ^ tr (r|L(p)rg_(p)t) = o , (27) 

which are straightforward generalizations of Eq. (pISl). Note that the operator vanishes 
in the Born approximation between quark states, which is another reason why we do not 
take it into account. 
Rewriting Eq. (^) as 

Ol\^.) = Z^Hai^) {o^Ha) + ^|^^^(«)) (28) 

we see that 0^^\jji) will have a multiplicative dependence on jj, (cf. Eq. (|T7|) ) only if the ratio 
Z'^ {ajj) / Z^^\ajj) does not depend on /x. The scale dependence will then completely reside 
in Zl^l. 



IV. SIMULATION DETAILS 

We have obtained numerical results for matrix elements and Z factors in quenched sim- 
ulations at /5 = 6/(7o = 6.0, 6.2, and 6.4 [g^ = bare coupling constant on the lattice). 
Whereas our original calculation [Q at /3 = 6.0 used Wilson fermions, we have meanwhile 
switched to non-perturbatively improved fermions (clover fermions) in order to reduce 0{a) 
effects. The value of the clover coefficient csw is taken from Ref . |^ . Since we have not 
improved the operators, there will still be residual 0{a) effects in the matrix elements and 
the renormalization factors. A few details of our computations are collected in Table |. 

The matrix elements are calculated from several hundred configurations for each /3. To 
compute the renormalization factors we use a momentum source ^ . Therefore the statistical 
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error is oc (^A^'conf)"^''^ for Nconi configurations on a lattice of volume V, and we already get 
small statistical uncertainties even from a small number of configurations, four in our case. 
(There is of course a price to be paid, the calculation for each momentum is independent, so 
the number of inversions of the fermion matrix is proportional to the number of momentum 
values.) The main source of statistical uncertainty in our final results is from the matrix 
elements, not the Z values. 

The momenta in the vertex functions used for the evaluation of the renormalization 
factors have been chosen close to the diagonal in the Brillouin zone in order to keep cut- 
off effects as small as possible. One should bear in mind that this diagonal extends up to 
= An'^/a^, but we use only momenta with < vr^/a^. 

In each case, the calculations are done at three (or more) values of the hopping param- 
eter K determining the bare quark mass so that we can extrapolate our results (both the 
bare matrix elements and renormalization factors) to the chiral limit. The extrapolation is 
performed linearly in m^, the square of the pion mass. 

The bare reduced matrix elements are calculated from three-point functions in the stan- 
dard fashion (see, e.g., Ref. p). In Eq.(25) of Ref. the ratios of three- to two-point 
functions for the 02 operator O^^^ and the d2 operator 0'^! are given. For the operators 

and (9°" the ratios and ratio factors are the same as for the d2 operator. The matrix 
elements are collated in Tables |T|, |ITT|, fV\ separately for u and d quarks in the proton. Here 
^2^^ and d2 correspond to the operators ( p2D and (|23|), respectively. In addition we give the 
pion masses (in lattice units) which we use in the chiral extrapolations. They are mostly 
taken from Ref. [^. Note that all our errors are purely statistical. They were determined 
by the jackknife procedure. 

V. NUMERICAL RESULTS FOR RENORMALIZATION COEFFICIENTS 

Let us begin the more detailed presentation of our numerical results with the renormal- 
ization factor Z^^^(afi) of the multiplicatively renormalizable 02 operator (|20|) . We convert 
our MOM numbers to the MS scheme using 1-loop continuum perturbation theory as de- 
scribed in P]. In Fig. |I] we show the fx dependence of Z^^^ extrapolated to the chiral limit. 
Results for Wilson fermions can be found in Ref. . Note that at scales /x^ exceeding a few 
times the lattice cut-off strong lattice artifacts may be present so that the corresponding 
results should not be taken too seriously. 

Turning to the more subtle renormalization of the ^2 operator (^2]) we must note that 
the conversion factor from our MOM scheme to the MS scheme has not yet been calcu- 
lated because of the complications caused by the mixing effects. Therefore we stick to 
the MOM numbers. Let us first consider the ratio Z'^{afi)/Z^^^{afi). As discussed above, 
Z'^ [ajj.) / Z^^^aji) should be independent of the renormalization scale n if the renormalized 
operator is to depend on /i multiplicatively. In Fig. |^ we show this ratio for our three j3 
values. It becomes approximately fiat for scales /i larger than about 3.5 GeV. While a scale 
of 3.5 GeV might seem to be somewhat too close to the cut-off for (3 = 6.0 and perhaps 
also for f3 = 6.2, it enters the region where lattice artifacts die out in the case of (3 = 6.4. 
Therefore we feel encouraged to apply the Z factors around /x = 3.5 GeV in order to evaluate 
structure function moments in the next section. 
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In Fig. ^we plot Z^^^{afi) in order to show the size of the multiphcative renormahzation 
in our approach. 



VI. NUMERICAL RESULTS FOR STRUCTURE FUNCTION MOMENTS 



Let us now discuss our nucleon matrix elements. In Figs. ^ we show the chiral extrap- 
olations of the bare values of a^^ and 0^2^' , respectively. Unfortunately, like in all other 
current QCD simulations our quark masses are rather large so that the extrapolation has to 
bridge quite some gap. A striking feature of the data is that the bare a'^^ values at /3 = 6.2 
are rather different from the values at the other two /9's. We can interpret this only as an 
unpleasantly large statistical fluctuation. The bare matrix elements in the chiral limit will 
be combined with the renormahzation factors of the preceding section to yield estimates of 
the renormalized matrix elements. 

We start with the twist-two matrix element 02. In the MS scheme with anticommuting 
75 the corresponding Wilson coefficient is given by (see, e.g., Ref. [p2| ) 



1 



X 



167r2 



7o/(2/3o) 



7l 

2/3o 



g\Q' 

167r2 



(29) 



with (3o = 11, pi = 102, 7o = 100/9, 71 = 141.78. (These are the numbers for Nf = flavors, 
appropriate for the quenched approximation.) The renormalized reduced matrix element is 
obtained from the bare value (extrapolated to the chiral limit) after multiplication with 
the non-perturbative renormahzation factor converted to the MS scheme. From Eq. (|^ 
we can then calculate dxx'^gi{x, Q"^). To avoid large logarithms in the Wilson coefficient 
we put = fi'^. In Fig. |^ we show the results for the proton and compare with the 
experimental value |^. While our lattice results at /3 = 6.0 agree surprisingly well with the 
experimental number, the above mentioned fluctuation makes them considerably larger at 
f3 = 6.2. Fortunately, they drop again at (3 = 6.4. For the neutron there is a similar effect, 
but due to the larger errors it is less significant. 

Let us now turn to our results for the twist-3 matrix elements. In this case it is unclear 
how to convert our MOM results to the MS scheme due to the mixing effects. Therefore 
we do not make use of the MS Wilson coefficient, which has recently been calculated ||2^ 
(with the 't Hooft-Veltman 75). It would change the final results for d2 by ^ 10%. Instead 
we use only the lowest-order approximation for the coefficient functions, i.e. the tree-level 
coefficients (|^), which are the same in all schemes, and define (by a slight abuse of notation) 



Ad) 



4") = + Q(")24'^) 



for the proton and the neutron, respectively. The renormalized values of ^2 / 
the proton are calculated from 



di^^ = z^'W^^^ + -z^d;^^K 



(30) 
(31) 

u, d in 



(32) 



8 



Remember that, besides the twist-three matrix element d2, dxx'^g2{x,Q'^) also contains 
a twist-two piece, the Wandzura-Wilczek contribution, see Eq. (||). To be consistent we 
restrict ourselves to the tree-level Wilson coefficients and the MOM matrix elements also in 
this contribution when computing dxx'^g2{x,Q'^) from (cf. Eqs. (|I]) and (^) 

/ dxx^5(2(x,(5^) = ^c/s - ^ / dxx^gi{x,Q^) . (33) 
Jo 6 3 Jo 

The moment dxx^g2{x,Q^) is plotted in Fig. |^ for the proton, where we have again 
identified = fi^. The experimental value is obtained by combining dxx'^gi{x, Q^) from 
Ref. [0 with d2 from Ref. [Q. Again we see the effect of the "fluctuation" at /? = 6.2. 

Comparing the proton results shown in Fig. ^ with the numbers presented in Fig. ^ one 
sees that dxx'^g2{x, Q^) is dominated by the twist-two operator. There is little room left 
for the twist-three operator, and one obtains rather small values for d2 as shown in Fig. ^ 
for the proton. In the neutron, 62 is even smaller in magnitude and hardly different from 
zero within the statistical errors. 

In Tables and ^ we present results at the scales 5 GeV^ and 10 GeV^, respectively. 



Note that here dxx'^gi{x,Q'^) includes the one-loop Wilson coefficient as well as the 
conversion factor to the MS scheme, both of which were neglected in the calculation of 
Jp dxx'^ g2{x , Q"^) for the reason explained above. The difference between the two sides of 
Eq. (|33D when evaluated with the numbers taken from the tables gives therefore an im- 
pression of the uncertainties originating from our incomplete knowledge of the perturbative 
corrections. 

In Figs. 1^, we fix the scale at 5 GeV^ and plot our results for the proton as well 
as for the neutron versus the lattice spacing a. Although an extrapolation to the continuum 
limit appears to be problematic, it is reassuring to see that we are getting close to the 
experimental numbers shown at a = 0. 

Of course, we should not forget that our computation suffers from various uncertainties. 
Apart from the fact that our treatment of the operator mixing is still incomplete, these 
concern, e.g., the influence of the quenched approximation, the extrapolation to the chiral 
limit, and the size of the lattice artifacts. Sea-quark effects are expected to be concentrated 
at small x, hence they should be suppressed by the factor x^ in the moment which we have 
considered. Therefore we may hope that the quenched approximation is reasonable in the 
case at hand. If indeed the valence quarks dominate, then it should also be justified to neglect 
flavor singlet contributions (like disconnected insertions and pure gluon operators), and it 
makes sense to consider proton and neutron matrix elements separately (as we have done) 
and not only flavor non-singlet combinations like ^2^'' — ^2"^ The quark mass dependence of 
our results is rather mild for the range of (relatively large) masses that we studied. Therefore 
the extrapolation to the chiral limit looks quite safe, although, of course, unexpectedly 
large effects at truly small masses cannot be excluded. Lattice artifacts are obvious in our 
renormalization factors (see, e.g.. Fig. |^). We have to expect them also in the nucleon matrix 
elements. Since we are working with non-perturbatively improved fermions we could in 
principle reduce their size by using improved operators. Unfortunately, the non-perturbative 
improvement of operators of the kind needed here is not straightforward and has yet to be 
worked out. In particular, improvement of the renormalization factors requires off-shell 
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improvement. Although there are some ideas on how to solve this non-trivial problem (see, 
e.g., Ref. 1 2^), an implementation for the operators considered here is beyond our present 
possibilities. 



VII. CONCLUSIONS 

In this paper we have tried to obtain a more reliable lattice estimate of the twist-three 
nucleon matrix element 62 improving on our first calculation in several respects. We have 
made a serious attempt to take into acccount the most important part of the operator mixing 
which occurs in this case, namely the mixing with lower-dimensional operators. This could 
only be done non-perturbatively and led to a significant change in the results for 62 moving 
them close to the experimental numbers. Thus the mixing with lower-dimensional operators 
seems to account for a large part of the difference between our previous computation and 
the experimental data. 

The calculations have been performed in the quenched approximation at three different 
values of 13 corresponding to three different values of the lattice spacing. While our results 
are still not good enough to allow for a meaningful extrapolation to the continuum limit, 
the mutual consistency of the values obtained for ^2 at the various /?'s indicates that dis- 
cretization effects are smaller than our statistical errors and corroborates our conclusion that 
the twist -three nucleon matrix element is rather small, in agreement with the experimental 
findings. 

We consider the computations at /5 = 6.4, i.e. at our smallest lattice spacing, to be most 
reliable. At = 5 GeV^ they yield the structure function moments (cf. Table |V]) 

X axx gi[^x) <^ Q QQQ ^ Q_QQ2 (neutron) 

for gi and 

Jo ^^^^> ~ \ -0.0002 ± 0.0017 (neutron) 
for g2. These numbers are to be compared with the experimental results [^,^ 

Cdxx'aAx) - I °-°^24 ± 0.0010 (proton) 

X cl^^ 9i[^) - I _o.0024 ± 0.0016 (neutron) ^"^^^ 

and 

Jo 92[x) - I Q_QQ29 ± 0.0035 (neutron) ^'^^^ 

respectively. 
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TABLES 



TABLE L Simulation parameters. In the third column ME indicates the calculation of nucleon 
matrix elements, whereas Z signifies the computation of renormalization factors. The lattice spacing 



a has been determined from the force scale tq |25] using tq = 0.5 fm; csw is the value of the clover 



coefficient. The matrix element calculations for the smallest quark mass {k = 0.1353) at /J = 6.4 
have been performed on a 323 X 54 lattice. 

/3 Lattice a^^[GeV] csw 

"eX) 16^ X 32 ME 2T2 1.769 

6.0 24? X 48 Z 2.12 1.769 

6.2 24^ X 48 ME 2.90 1.614 

6.2 24^ Z 2.90 1.614 

6.4 32^ X 48 ME 3.85 1.526 

6.4 32^ X 40 Z 3.85 1.526 



TABLE IL The unrenormalized reduced matrix elements a2, ^2 ^-iid for u and d quarks in 



the proton 


at /3 = 6.0. 


Also given are 


the pion masses. 










0.132 


0.1324 


K 

0.1333 


0.1338 


0.1342 




amir 


0.5412(9) 


0.5042(7) 


0.4122(9) 


0.3549(12) 


0.3012(10) 


0.0 




0.114(8) 
-0.029(3) 


0.114(8) 
-0.032(4) 


0.107(11) 
-0.037(7) 


0.09(2) 
-0.032(12) 


0.08(2) 
-0.047(16) 


0.08(2) 
-0.046(11) 


^[5](d) 


0.0063(12) 
-0.0041(6) 


0.0028(13) 
-0.0027(7) 


-0.010(2) 
-0.0008(14) 


-0.023(5) 
-0.005(2) 


-0.028(6) 
-0.001(4) 


-0.037(4) 
0.001(2) 




-0.216(12) 
0.050(4) 


-0.228(12) 
0.050(3) 


-0.246(18) 
0.052(6) 


-0.29(4) 
0.064(12) 


-0.27(3) 
0.044(17) 


-0.30(3) 
0.057(12) 



TABLE III. The unrenormalized reduced matrix elements 02, ^2 ^^i^l d!^ for u and d quarks 



in the proton 


at P = 6.2. Also 


given are the 


pion masses. 








0.1333 


0.1339 


0.1344 


0.1349 






0.4136(6) 


0.3570(10) 


0.3034(6) 


0.2431(7) 


0.0 


4"^ 


0.142(10) 
-0.033(4) 


0.137(15) 
-0.030(5) 


0.157(17) 
-0.034(9) 


0.16(3) 
-0.031(15) 


0.17(3) 
-0.030(13) 


45] (.) 


-0.0017(14) 
-0.0027(6) 


-0.017(2) 
0.0006(11) 


-0.031(5) 
-0.0003(17) 


-0.051(11) 
-0.000(4) 


-0.065(7) 
0.004(3) 




-0.34(2) 
0.067(6) 


-0.37(3) 
0.065(8) 


-0.44(4) 
0.072(13) 


-0.50(7) 
0.07(3) 


-0.54(7) 
0.07(2) 
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TABLE IV. The unrenormalized reduced matrix elements 02 , 02 and for and d quarks in 
the proton at /3 = 6.4. Also given are the pion masses. 



(1 "1 ■J s; 




n 1 'J 1 ( ; 

U. i^) l^> 


U , 1 o ■. ) 






arrin 0.3213(8) 


0.2836(9) 


0.2402(8) 


0.1933(7) 


0.1507(8) 


0.0 


o^"' 0.123(8) 
a^^ -0.030(4) 


0.092(13) 
-0.018(5) 


0.114(13) 
-0.032(7) 


0.102(17) 
-0.033(10) 


0.14(4) 
-0.01(2) 


0.096(19) 
-0.023(10) 


4^^' '"^ -0.0149(14) 
4^' '''^ 0.0014(7) 


-0.024(3) 
0.0021(13) 


-0.038(4) 
0.0062(17) 


-0.055(6) 
0.012(3) 


-0.055(15) 
0.015(8) 


-0.068(6) 
0.013(3) 


dj'^Va -0.40(2) 
4'-^^ /a 0.088(6) 


-0.41(4) 
0.075(12) 


-0.47(4) 
0.105(13) 


-0.51(5) 
0.12(2) 


-0.55(13) 
0.15(5) 


-0.56(5) 
0.12(2) 




TABLE V. Results for /x^ = 


(5^ = 5GeV^. 








6.0 




6.2 




6.4 


jj^ dxx'^g^^ (x) 
!idxx^g^\x) 


0.046(13) 
0.008(4) 
0.012(3) 
-0.007(3) 




0.11(2) 
0.017(9) 
0.029(5) 
-0.019(4) 




0.066(14) 
0.017(7) 
0.017(4) 
-0.010(3) 


/Jda;a;2yJ'^(x) 


-0.017(8) 
-0.003(2) 
-0.004(2) 
0.0026(16) 




0.009(11) 
-0.001(4) 

0.002(3) 
-0.002(2) 




0.000(8) 
-0.001(3) 

0.000(2) 
-0.0002(17) 




TABLE VL Results for /x^ = 


= lOGeV. 








6.0 




P 

6.2 




6.4 


4^) 

jj^ dajx^g-j^^ (a;) 
^ldxx^g^^{x) 


0.040(12) 
0.002(4) 
0.010(3) 
-0.007(2) 




0.098(18) 
0.006(7) 
0.025(5) 
-0.017(4) 




0.057(12) 
0.008(6) 
0.015(3) 
-0.009(3) 


4"^ 

/o^dxa:;25i")(x) 
/(Jdxx24"^(a;) 


-0.015(7) 
-0.003(2) 
-0.0037(18) 
0.0021(14) 




0.007(9) 
-0.002(3) 

0.002(2) 
-0.0016(18) 




0.000(7) 
-0.001(3) 

0.0001(19) 
-0.0002(14) 
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FIG. 1. The renormalization constant Z^^^ in the MS scheme. 
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FIG. 2. The ratio Z'^(a//)/Z[5] (a/x). 
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FIG. 3. The renormalization factor Z^^^a/j,). 
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FIG. 4. Chiral extrapolation of the bare matrix element in the proton 
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FIG. 5. Chiral extrapolation of the bare matrix element ^2 in the proton. 
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FIG. 6. The moment Jq dxx^gi{x,Q'^) for the proton. The square indicates the experimental 
value S. 



20 



0.01 



0.00 



-0.01 - 



-0.02 



-0.03 



0.00 



-0.01 



-0.02 



-0.03 - 



0.00 



-0.01 



-0.02 



-0.03 



-0.04 



/3= 6.40 




4 5 6 7 

Q (GeV) 



10 



FIG. 7. The moment Jq dxx^g2{x,Q'^) for the proton. The square indicates the experimental 
value obtained by combining results from Q and Q. 
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FIG. 9. The moment Jq dxx^gi{x, Q"^) at = 5 GeV'^ for the proton (open symbols) and the 
neutron (filled symbols) plotted versus the lattice spacing a. The squares at a = indicate the 
experimental values [^. 
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FIG. 10. The moment Jq dxx'^g2{x, Q"^) at = 5 GeV^ for the proton (open symbols) and the 
neutron (filled symbols) plotted versus the lattice spacing a. The squares at a = indicate the 
experimental values obtained by combining results from and [Q]. 
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FIG. 11. The reduced matrix element d2 at fx^ = 5 GeV^ for the proton (open symbols) and 
the neutron (filled symbols) plotted versus the lattice spacing a. The squares at a = indicate the 
experimental values Q. They are plotted with a slight horizontal offset to avoid overlapping error 
bars. 
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